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Abstract 

Using the facilities of the KASCADE Central Detector EAS muon arrival time 
distributions, observed with reference to the arrival time of the first locally registered 
muon, and their correlations with other EAS observables have been experimentally 
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investigated. The variation of adequately denned time parameters with the dis- 
tance from the EAS axis has been measured. The experimental data enable a 
study of the sensitivity of such local arrival time distributions, which characterise 
the structure of the shower disc, to the mass composition of cosmic rays in the 
energy region around the knee. For that purpose, nonparametric multivariate even- 
by-event analyses have been performed for an estimate of the mass composition 
specified by three different mass groups, invoking detailed Monte Carlo simulations 
of the EAS development. It turns out that local muon arrival time distributions, 
without information on the curvature of the shower disc, display a minor sensitivity 
to the mass of the EAS inducing particle, at least for distances from the shower axis 
< 100 m. The measurements comprise a subset of all EAS events registered by 
K AS CADE due to the observation conditions of the arrival time distributions, with 
a threshold of the muon energy Efh = 2.4 GeV and a minimum multiplicity nth for 
being accepted in the observed data samples. This subset is sensitive to variations 
of the integral EAS muon energy spectrum. By studying the event acceptance in 
the registered samples on basis of Monte Carlo simulations a test of the consistency 
of the Monte Carlo simulations with the data is enabled, comparing the results in- 
ferred from observations at different R^ and different nth values. Within the present 
uncertainties the results of such a test show a remarkable agreement of the expe- 
rimental findings with the Monte Carlo simulations, using the QGSJET model as 
generator of the high-energy hadronic interactions. 

Key words: Extensive Air Showers, Arrival Time Distributions, Primary Mass 
Composition, Model Test 
PACS: 98.70Sa, 96.40Pq 



1 Introduction 

The time delay of the particles inside the front of Extensive Air Showers 
(EAS) and the temporal structure of different EAS components are a subject 
of longstanding interest of cosmic ray research. In fact the first experimen- 
tal studies have been performed in 1953 by Bassi, Clark and Rossi [1] and 
Jelly and Whitehouse [2], followed by many others [3]-[ll]. The renewed in- 
terest arises from recent measurements using advanced detector facilities like 
the COVER PLASTEX detector within the GREX array [12] -[14], the EAS 
TOP array [15] or the facilities of the KASCADE experiment [16]. Investi- 
gations of the structure of the EAS disc by the KASCADE experiment are 
focused in particular on the muon component [17]-[21]. Due to the reduced 
influence of multiple Coulomb scattering of GeV muons and the absence of 
absorption, muons travel nearly in straight lines from the locus of production 
to the observing detector. Thus muon arrival time distributions, observed at 
large distances from the shower axis, are expected to map the longitudinal 
EAS development and to reflect the distribution of production heights via the 
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time-of-flight of the muons. With some simplifications, using a triangulation 
procedure, the distribution of the production heights could be estimated from 
the time delay of the observed muons relative to the arrival time of the shower 
centre [17,22,23]. Basically the same information could be alternatively de- 
duced from the distribution of angles of muon incidence relative to the shower 
axis [24]. Corresponding measurements are also a current subject of KAS- 
CADE investigations using the Muon Tracking Device [25]. Angle-of-incidence 
observations or the combination with arrival time measurements (as discussed 
with the term Time Track Complementary [22,26]) do not reveal additional 
basic information [23] though they provide an interesting practical alternative. 

The basic sensitivity of muon arrival time distributions (whose phenomeno- 
logical features seen with the KASCADE experiment are reported in Ref. [20]) 
to longitudinal EAS development and to the elongation rate [27-29] puts the 
question which particular parameters of the observed distributions provide 
useful signatures of the mass of the EAS primary and under which conditions 
they are helpful for the determination of the mass composition of the primary 
cosmic rays. The present paper addresses experimentally the question of the 
sensitivity to the primary mass on the basis of results from the KASCADE ex- 
periment, in particular of EAS event-by-event measurements of the temporal 
EAS structure (shower thickness) observed at relatively small distances to the 
shower axis < 100 m for the primary energy range around the knee. Me- 
thodically non-parametric statistical analysis techniques of multivariate dis- 
tributions [30] are applied for the mass classification of the observed EAS by 
their time parameters, derived from the muon arrival time distributions, and 
by their correlations with other EAS observables (like the shower size N e and 
the muon content N^). For muon arrival time distributions the methods are 
outlined in Refs. [17,23]. These procedures avoid the bias of pre-chosen func- 
tional forms and compare to reference patterns, derived from extensive Monte 
Carlo simulations, by use of the EAS Monte Carlo simulation code CORSIKA 
[31], thus including all natural EAS fluctuations. By applying additional de- 
tailed detector simulations based on GEANT [32], the distortions arising from 
the experimental conditions [33] are taken into account. 
The analysis of the EAS observations in terms of the mass of the primary intro- 
duces necessarily some model dependence by the high-energy hadronic inter- 
action models invoked as generators of the Monte Carlo simulations. Therefore 
the results are subject of the uncertainities of the hadronic interaction model, 
used for the analysis. A possible way to evaluate the quality of a model is to 
derive the primary mass composition from the analysis of different observables 
and to consider the agreement or disagreement of the results. This idea has 
been worked out in a consistent and efficient manner [34] with an extended 
set of observables of the KASCADE experiment, showing that there exist 
systematic differences in the estimate of (InA), e.g., derived from different 
combinations of (correlated) observables [35], revealing the limitations of po- 
pular models like QGSJET [36], VENUS [37] and SIBYLL [38]. This indicates 
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that improvements of the models are urgently needed on basis of experimental 
findings. The QGSJET model is actually considered to be that one of the best 
internal consistency [39]. Hence we refer furtheron to the QGSJET model. 
The muon arrival time measurements imply a particular selection of EAS, 
with a distortion of the original mass composition in the subset of classified 
EAS. This is due to the energy threshold E th (2.4 GeV for KASCADE) of 
the muon detection and the condition of a minimum multiplicity nth (= 3 
muons detected by the KASCADE Central Detector) for reconstructing an 
arrival time distribution. The distortion is dependent on the particular values 
of E th , n t h, on the mass A of the primary and the muon content N^, but 
also on the radial distance of the arrival time measuring detector from the 
shower axis. In order to restore the original mass composition corresponding 
efficiency (acceptance) correction factors have to be applied to the identified 
mass groups. These factors can be only determined by Monte Carlo simula- 
tions, studying the cuts implied by the experiment and their effect on the 
efficiency. An efficient and sensitive test of the interaction model used for the 
simulations and of the particle tracking procedures is enabled by applying the 
calculated correction factors to the results found at different and to look 
for the consistency of the resulting mass composition. Such a test of the Monte 
Carlo simulations and of the QGSJET model, especially in view of predictions 
of the muon energy spectra and muon densities and their fluctuations in event- 
by-event observations, is an important aspect of muon arrival time studies of 
this paper. The results add to the conclusions about the difficulties of the 
model in interpreting consistently muon density measurements with different 
muon energy thresholds [40]. 



2 Muon arrival time measurements 

Measured arrival times < t 2m < t 3m < ... of muons registered by timing 
detectors at a distance from the shower axis, corrected by ±-R M tan#/c 
(c - speed of light) in order to eliminate the distortions due to the shower 
inclination, refer to an experimentally provided zero-time. Depending on the 
choice of the kind of zero-time reference, there are two different types of muon 
arrival time distributions. By the use of the arrival time r c of the shower core 
as reference global arrival times are observed [33]: 



(1) 



e.g. for the foremost muon registered at R^. 



(2) 
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This type of time distributions informs about the curvature of the shower 
front as well as about the structure of the disc. The arrival time of the shower 
core is difficult to determine with sufficient precision. Therefore the analysis 
has been restricted to " local" times, which refer to the arrival of the foremost 
muon locally registered by the detector: 



Art *W = n»(R») - r^R,) . (3) 

(with omitting further the label loc), informing only about the thickness and 
the structure of the muon disc. Implications of observations of local muon 
arrival time distributions due to the fluctuations of the arrival of the first 
registered muon have been discussed in Ref . [33] . For event-by-event observa- 
tions with a fluctuating number of muons (multiplicity), the individual relative 
arrival time distributions can be characterised by the mean values Ar mean , and 
by various quantiles Ar q , like the median Aro.50, the first quartile Aro.25 and 
the third quartile Aro.75 (for details of the definition see Refs. [19-21]). Their 
mean values and dispersion (standard deviations) represent the time profile of 
the EAS muon component. 



KASCADE, whose layout is described in more detail in Ref. [16], is a mul- 
tidetector system, installed in Forschungszentrum Karlsruhe (110 m a.s.L), 
Germany, for the observation of extensive air showers in the primary energy 
range around the knee. One part is an array of 252 detector stations, dis- 
tributed over an area of 200 x 200 m 2 on a grid of 13 m spacing for measuring 
the electron-photon component and the muon component with a threshold of 
5 MeV and 230 MeV, respectively and providing the basic information about 
arrival direction, core position, electron and muon sizes of the observed EAS. 
In particular, from the data of the field array the so-called truncated muon 
number iV* r , i.e. the muon density integrated between 40 m and 200 m, is 
derived and used in the KASCADE case (between 10 14 and 10 16 eV) as an ap- 
proximate mass independent energy identifier [41]. The location of the shower 
core can be determined (inside the fiducial area) with an uncertainty of less 
than 3 m. The arrival direction of the shower is reconstructed with an uncer- 
tainty better than 0.5°. Details of the reconstruction procedures are described 
in Ref. [42]. As additional muon detector Limited Streamer Tubes tracking 
detectors have been installed in an underground tunnel for measurements of 
the lateral distribution (E t h = 0.8 GeV) and of muon angles-of-incidence dis- 
tributions [25]. 

The muon arrival time measurements use, in particular, the facilities of the 
Central Detector [43] of KASCADE. It is basically an iron sampling calorime- 
ter (with an area of 16 x 20 m 2 , set up with liquid TMS and TMP ionisation 
chambers) for the identification and energy measurement of hadrons. In the 
basement of the set-up, below 3800 t of iron and concrete, large-area position- 
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sensitive multiwire proportional chambers (MWPC) [44] are operated for the 
identification of muons with 2.4 GeV energy threshold. The performance of 
the detection system of the MWPC is improved by a layer of streamer tubes 
[45]. The trigger plane of the calorimeter is a system of 456 plastic scintil- 
lation detector elements (47.5 x 47.5 x 3 cm 3 in size, each separated by a 
wavelength shifter) for providing a fast trigger signal (in addition to the trig- 
ger from the field array) for the MWPC, and for the timing measurements. 
Fast electronics [46] records low (muons) and high energy deposits (cascading 
hadrons). In order to remove the signals from cascading hadrons an upper 
limit of the energy deposit of 30 MeV in each of the 456 scintillation coun- 
ters is imposed. An amount of 24 millions of KASCADE events are used for 
the analysis [47]. Muon arrival time distributions have been reconstructed for 
muons with E th = 2.4 GeV and for events with a number n > 3 of registe- 
red muons; after applying also some general cuts concerning the core position 
(within 100 m from the array centre), the angle of EAS incidence (< 40°) 
and logioNjJ' (> 3.6) the sample shrinked to approx. 240 000 showers. The 
phenomenological features of the observed muon arrival time distributions 
have already been reported in Ref. [21], where further experimental details 
are communicated. 



3 EAS simulations 



The interpretation of the measured muon arrival time distributions and their 
correlations with other shower parameters need a-priori knowledge to be de- 
duced from Monte Carlo simulations of the EAS development. The present 
analysis is based on simulations with the code CORSIKA (version 5.62) with 
a full and detailed simulation of the detector response. The simulations use the 
QGSJET (version 1998) model [36] as generator for high-energy interactions 
and GHEISHA [48] for interactions below Ei a b = 80 GeV. The electromagnetic 
part is treated by the EGS4 program [49]. Earth magnetic field, observation 
level and particle detection thresholds have been chosen in accordance to the 
experimental situation. The U.S. standard atmosphere [31] has been adopted 
for the simulations. The simulations have been performed for three different 
classes of primaries: protons (H) for the light group, oxygen (O) for the CNO- 
group and iron (Fe) for the heavy group. The energy range covered by the 
simulations extends from 5.0 • 10 14 eV to 3.06 • 10 16 eV and the zenith angles 
comprise the range of 0° < 9 < 40°. The centres of the simulated showers 
have been positioned inside a quadratic area of 210 x 210 m 2 , slightly extend- 
ing the KASCADE area. For each primary type approximately 90000 showers 
have been simulated, with decreasing numbers of shower events for higher 
primary energies and the larger zenith angles, due to restrictions in the com- 
puting time. Finally a sample for 0° < 9 < 24° (corresponding to a third of 
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Fig. 1. The distribution of various shower observables of the sample, prepared by 
simulations for different types of primaries, compared with the distribution of the 
sample of experimental data for 4.05 < logioNf[ < 4.28, 80 m < Rn < 90 m and 



0° < < 24°. The quantity n is the multiplicity of muons (Eth 
with the MWPC setup of the Central Detector. 



=2.4 GeV) detected 



the considered sec 9 - range of (1, sec 40°)) was used in the analysis. The simu- 
lated primaries (separately for each of the cases H, O, Fe) have been weighted 
according to: 



dN = const. ■ E 7 • sin 9 ■ cos 9 ■ D core ■ dEo ■ d9 ■ dD core (4) 

where dN is the number of primaries with energy between E , E + dE , zenith 
angle between 9,9 + d9, intersecting the plane of the detector array between 
D core , D core + dD core . A fixed spectral index 7 = 2.7 over the whole primary 
energy range is adopted for all primaries (H, O, Fe). More details about the 
simulations are given in Refs. [50] and [47]. As input for the further analysis 
of the muon arrival time distributions and their comparison with model pre- 
dictions the following shower parameters and their correlations are regarded: 



7 



• the shower size N e ; 

• the truncated muon number iV* r ; 

• the multiplicity n of muons (E th = 2.4 GeV) detected in the facilities of the 
Central Detector; 

• the quantiles Ar q of the local muon (E th = 2.4 GeV) arrival time distribu- 
tions at various distances from the shower centre; 

• the reduced quantiles At* = Ar g / p M of the local muon arrival time distri- 
butions (E th = 2.4 GeV), i.e. the quantiles divided by the density p M , where 
p M is estimated from the observed multiplicity and the effective area of the 
muon detector set-up. 

Fig. 1 displays frequency distributions of some shower parameters, prepared 
by the simulation for different primaries and compared to the actual experi- 
mental observations for particular logwNj[ and bins. Each distribution is 
normalised (only for the presentation of N e the distribution of oxygen- induced 
EAS is additionally shown). There are some discrepancies in the multiplicity 
distribution between simulations and experimental observations. They may 
arise from an imperfect adoption of the primary energy spectrum and from ef- 
ficiency effects not fully corrected for. Since they also affect the fluctuations of 
the local time parameters originating from the multiplicity dependence (dis- 
cussed in detail in Ref. [33]), this analysis uses preferentially the reduced 
quantities At*. The fluctuations largely cancel in the reduced parameters. 
Actually it has been argued that the reduced quantiles At* could exhibit en- 
hanced mass sensitivity since they include approximately the A.r q {R^)—p^{R^) 
correlation [51]. However this sensitivity is obscured by the limitations of the 
experimental response. CORSIKA simulations show [23,51] that the age pa- 
rameter is carrying relevant information about the longitudinal development. 
This aspect has not been worked out by the present analysis. 



4 Non-parametric multivariate analyses 

Non-parametric statistical methods are applied in studies of multidimensional 
distributions of observables allocating the single observed events to different 
classes (in our case to proton, oxygen, or iron primaries) by comparing the 
observed events with the simulated distributions without using a pre-chosen 
parameterisation. The methods of decision making and procedures of appli- 
cations to cosmic ray data analyses are extensively described in Refs. [30,35] 
and outlined for analyses of muon arrival time distributions in Refs. [17,23]. 
The procedures take into account the effects of the EAS fluctuations in a 
quite natural way and are able to specify the uncertainties, by an estimate 
of the true-classification and misclassification probabilities. The classification 
probabilities are determined by the extent to which the likelihood functions 
of the individual classes, derived from the simulations, are overlapping. Ba- 
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sically such pattern recognition methods, using trained neural networks or 
Bayes decision rules, depend on the hadronic model generating the reference 
patterns for the experimental observables. The concept of the present analysis 
is to determine in a first step the sensitivity of various observable correlations 
to the mass of the primary, by applying the one-leave-out test to the sample 
of distributions prepared by Monte Carlo simulations (QGSJET model). The 
one- leave-out test determines the probability that a multidimensional event, 
taken from the considered (simulated) distribution, will be correctly ("true") 
or incorrectly ("false") classified by the pattern recognition procedure (see ap- 
pendix of Ref. [35]). With respect to the experimental data sample the studies 
are performed not only for six different logioNj[ ranges between 3.60 and 
5.00 (covering the knee-range) but also for four different .R^-ranges (45-65 m, 
65-80 m, 80-90 m, 90-100 m). This allows to explore trends of the classifi- 
cation and misclassification probabilities and other quantities with varying 
Using the determined classification probabilities the mass composition is 
reconstructed from the data samples measured at different ^-ranges. The 
resulting mass composition generally differs from the primary mass composi- 
tion searched for, since the analysed data samples represent selections in the 
multidimensional space of all EAS observables. Hence efficiency correction fac- 
tors have to be determined, whose application should lead to a primary mass 
composition consistent for all -R^-ranges. The determination of the efficiency 
correction factors invokes again Monte Carlo simulations, and the extent of 
agreement of the primary mass composition resulting from different ranges 
provides a test for the procedures and the generator of the Monte Carlo simu- 
lations. 



4-1 The classification and misclassification probabilities 

As example in Fig. 2 the true-classification and misclassification probabilities, 
inferred with the one-leave-out test from the simulation sample, are displayed 
for H, O, Fe primaries and for different combinations of observables of EAS 
events. The events have been registered with the conditions of muon arrival 
time measurements in a particular {log w Nj[, R^} - range. With little surprise 
it is noted that the mass discrimination is dominated by the {NjJ', N e } corre- 
lation (see Fig. 3), and it is also evident that local time parameters (shown 
for the median Aro.50 or the third quartile Aro.75 and their reduced values 
Atq Atq 75 ) have minor influence on the discrimination. We emphasise that 
this statement holds for local arrival time distributions in the studied range of 
relatively small distances R^ (< 100 m) from the shower centre. Theoretical 
studies [47,51] indicate more pronounced differences in the temporal structure 
of the shower disc with increasing R^ and logioNf[. In contrast to results ig- 
noring the detector response [51] the use of reduced parameters Atq 50 , Ato* 75 
does not significantly improve the mass discrimination. Fig. 3 shows the true 
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Fig. 2. Classification and misclassification probabilities determined by Bayes deci- 
sion making for three different classes (H, O, Fe) and for distribution of different 
types of correlated EAS observables with 3.83 < log±oNj[ < 4.05, 80 m < < 
90 m, 0° < 8 < 24°. The results of different observable combinations are connected 
by (full and dashed) lines. 
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Fig. 3. The dependence of the averaged true-classification probability from {Nj[ , N e } 
and {N* r ,N e , Ar* 50 , Ar* 75 } correlations (0° < 6 < 24°). 
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classification probabilities (and uncertainties) averaged over all classes with in- 
creasing logioNj[ indicating a marginal systematic improvement if the time in- 
formation is added. The uncertainties are determined by the bootstrap method 
[30] , which basically consists of applying the classification procedure for a test 
sample several times, thus deriving an average value and the variance. 

4-2 Reconstruction of the mass composition from the observed data samples 

The true-classification P^i and misclassification -Pj_>j probabilities, deduced 
for all studied {logioNff, R^} ranges, are used for the reconstruction [23] of the 
mass composition of the samples of registered events by inverting the system 
of linear equations: 

n' H = n H • Ph->h + n • Po->h + n Fe ■ P Fe ^ H 

n' = n H ■ P H ->o + n • Po-^o + n Fe ■ P Fe ^o ( 5 ) 

n Fe = n H ■ PH^Fe + n • Po~*Fe + n Fe • P Fe ^Fe 

where nn, no, n Fe are the true numbers, defining the mass composition in 
the sample N = nu + no +n Fe , getting altered to n' H , n' Q , n' Fe because of the 
misclassifications. In Fig. 4 the reconstruction of an arbitrarily adopted mass 




3.8 4 4.2 4.4 4.6 

Fig. 4. Test of the reconstruction of an arbitrarily adopted mass composition of the 
event samples on basis of the {Nff , N e , Atq 50 , Atq 75 } correlation, observed at 80 m 
< R„ < 90 m, 0° < 9 < 24°. 
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composition of the event samples (displayed by the corresponding symbols), 
is shown with the resulting uncertainties as varying with logioNj[ . The result 
of the application of the reconstruction procedures to experimental samples 
measured with KASCADE, observed at different distances from the shower 
centre, are shown in the upper panels of Fig. 5. The uncertainties are influenced 
by the limited statistical accuracy of the data samples, in particular for the 
bins with larger logioNff values i.e. primary energies. Nevertheless it is obvious 
that mass compositions of measured KASCADE samples differ for different 
bins, since the observation conditions lead to mass dependent differences in 
the observation efficiency at different R^, leading to distortions of the deduced 
primary mass composition. This feature should be removed after applying a 
correct efficiency correction procedure. That necessarily implies again the use 
of the particular hadronic interaction model. 



4-3 Reconstruction of the primary mass composition 



Efficiency correction factors (Ch, Co, CVe) have been calculated in order to 
adjust the mass composition of measured KASCADE samples (Ph, Po, PFe) 
to the primary mass composition (Pr, Pq, -P_Fe) according to the relation [47]: 



p* . p* . p* _ P H . Po . PFe / r n 

^O <^Fe 

As a first step, the simulated spectra given by eq. 4 have been normalised to 
the same value (=1000) for each type of primaries A (H, O, Fe), and for all 
simulated cases of E , 6 and core position ranges. For a given primary A, the 
detected spectra will appear distorted at the ground level due to the absorption 
in the atmosphere and the selection cuts. Only a fraction of the original events 
are detected by the KASCADE detector and reconstructed successfully. The 
major influence on the values of the efficiency correction factors arises from 
the applied cuts on the reconstructed shower events. The correction factors 
depend on the R^, 9 and log w N^ ranges, on the multiplicity threshold for the 
2.4 GeV muons and, of course, on the primary type. For a given primary A, in 
a certain R^+AR^ 6,6+ AO; log 10 N* ,log w Nf[ +Alog 10 Nj[; n,n+An} 
(multidimensional) bin, the correction factor {C a) is given by the sum of the 
weights (wi) of the p simulated events which are accepted showers in the 
multidimensional bin: 



CU = (7) 
1=1 

The weights are determined via the normalised mass spectra (eq. 4) which 
get filtered by the observation conditions. Considering the equivalent number 
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Fig. 5. Variation of the mass composition of measured KASCADE samples with 
logioNjf estimated from the {Nff, N e , Atq 50 , Atq 75 } correlation and the primary 
mass composition i.e. before (samp, comp.) and after (corr. comp.) efficiency cor- 
rection shown for different ranges from the shower centre, for the zenith angle 
range 0° < 6 < 24°. 
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Fig. 6. Variation of the efficiency corrected mass composition and of (InA) with 
for various logioNjf ranges; 0° < < 24°. 



of events q [52] given by: 
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•t=i 

the statistical uncertainty (error of the mean value) has been calculated for 
C A by: 



5C A = C A -^- = ^. (9) 

Fig. 5 displays for different ^-ranges the /ogio-/V"* r -variation of the mass com- 
position of the measured KASCADE samples extracted from the {Nj[,N e , 
Atq 50 ,Atq 75 } correlation and the primary mass composition resulting after 
the correction concerning the biased acceptance by the specific observational 
conditions. In the lower part of Fig. 6 the results are shown as variation of 
(InA) with for different log 10 Nj[ ranges. The uncertainties of the correc- 



45m<R^65m 65m<R^80m 80m<R„^90m 90m <R„< 100m 
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Fig. 7. Variation of the calculated acceptance correction factors with for four 
i? M ranges and three log±oNf[ ranges in the zenithal angle range of 0° < 9 < 24°. 
The factor S to t/S r i ng (S tot = 210 x 210 m 2 ) accounts for the different geometrical 
areas for different ranges. 
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tion factors (eq. 9) have been included by propagation via eq. 6. It should 
be stressed that the displayed error bands do not include the uncertainty of 
the adopted spectral indices 7. Within the overall uncertainties the results 
obtained for different ranges are in fair agreement. This can be conside- 
red as a proof of the consistency of the Monte Carlo simulations invoked for 
the analysis of the data. The correction factors Ca give the fraction of the 
primary energy spectrum of the considered primary mass which contributes 
to the {Nj[,N e , ...} sample observed with conditions specified by the energy 
threshold of the registered muons E th , the distance from the shower centre 
and the multiplicity threshold n t h- In order to give an impression on these fac- 
tors, they are shown for various cases with their variation with n t h in Fig. 7. 
The variation of Ca is more pronounced at smaller Nj[ values (i.e. primary 
energies), and as expected from the muon lateral distributions the acceptance 
of Fe events in the observation sample is higher than for proton events. 



5 Test of the consistency of Monte Carlo simulations 

As emphasised above, the event samples registered for these measurements 
are particular selections of all EAS events, determined by the observation 
conditions, especially by the energy threshold of the detected muons and the 
multiplicity threshold. Thus the event selections are affected by the lateral 
variation and mass dependence of the muon energy spectrum of EAS. The 
mass composition of these samples, as inferred from the correlations of the 
observed EAS observables, are altered as compared to the primary composi- 
tion, in a way varying with the distance from shower centre. The acceptance 
or efficiency factors Ca for correcting this effect depend on the Monte Carlo 
simulations and their ingredients like the model descriptions of the hadronic 
interaction. Since the measurements are performed at different distances from 
the EAS centre, the variation of the results after applying the Ca factors im- 
plies a consistency test of the Monte Carlo simulations since the final result 
should be independent from within the given uncertainties. This test can 
be refined and be more stringent when not only the variation with R^, but also 
the variation with the multiplicity threshold n t h is scrutinised. In Fig. 8 the 
results of such a refined test are shown, based on the dominant {iV* r , N e } cor- 
relation within the observed samples. Observables describing the (local) EAS 
time structure (which are shown to be of minor sensitivity) have not been 
included in the classification procedure though their measurement defines the 
test sample. The results have been obtained by repeating the full data ana- 
lysis, described in the previous sections, for different multiplicity thresholds. 
In view of the rather specific observation conditions, sensitive to variations 
of the EAS muon energy spectrum, we consider the results, even if display- 
ing larger uncertainties, as a remarkable confirmation of the consistency of 
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Fig. 8. Variation of the deduced primary mass composition with the multiplicity 
threshold nth using the {Nj[,N e } correlation for the non-parametric classification 
procedure of the events observed at different distances from the shower axis. Shown 
are results for two different logxoNff ranges (corresponding to primary energies 
below and above the knee in the primary cosmic ray spectrum). 

the Monte Carlo simulations performed with the CORSIKA code using the 
QGSJET model. The test could be improved when samples resulting from spe- 
cific cuts would be analysed by simultaneously classifying the primary energy 
and the mass, using efficient pattern recognition methods as applied in other 
KASCADE studies [34,35]. The mass composition resulting from the present 
studies is in fair agreement with the results of those investigations. In parti- 
cular it corroborates the finding of an increase of {InA) beyond the knee, i.e. 
for log 10 Nj[ > ca. 4.2. 
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6 Conclusions 



The present investigation considers experimental data of the KASCADE ex- 
periment which are based on observations of local muon arrival time distri- 
butions, representing the variation of the EAS disc thickness with the dis- 
tance from the shower axis. The first focus of the studies is the question to 
which extent these time quantities and their correlations with other EAS ob- 
servables can help to improve the mass discrimination in the event samples 
observed with realistic experimental, i.e. with KASCADE conditions. Metho- 
dically the event-by-event analyses apply a non-parametric approach with ref- 
erence patterns provided by CORSIKA simulations of the EAS development, 
using as generators the QGSJET model for the high energy interaction part 
and GHEISHA for the low energy part. It turns out that similarly to other EAS 
observables, the local time quantities provide only a marginal contribution to 
the mass discrimination as compared to the dominant {Nff, N e } correlation, 
at least in the relatively limited range of < 100 m. It should be emphasised 
that this statement does by no means disqualify the sensitivity of global muon 
arrival time distributions measured relative to the arrival time of the shower 
centre. In fact, studies including the curvature of the shower front indicate a 
considerable enhancement of the mass discrimination, in particular at larger 
distances from the shower core and higher primary energies [47,51]. 
The event sample collected with the observation of muon arrival time distri- 
butions in the KASCADE experiment is a subset of all EAS events measured 
with 100% efficiency for log 10 Nf[ > approx. 3.6. As consequence of the muon 
energy threshold E t h = 2.4 GeV and the condition of a registered multipli- 
city n > 3 the original primary mass composition in the registered samples 
gets distorted in a way which is sensitive to lateral variations of the integral 
EAS muon energy spectrum. In order to determine the fraction of all shower 
events which are accepted in the specific subset, specified by the measuring 
conditions, Monte Carlo simulations have to be invoked for the calculation of 
corresponding acceptance or efficiency factors. The variation of such correc- 
tions with the distance of observation i? M from the shower axis and with the 
multiplicity threshold n t h provides the possibility to test the consistency of the 
Monte Carlo simulations with the data. Even admitting the large uncertain- 
ties due to the limited number of events observed and simulated, our results 
indicate a remarkable consistency of the performed Monte Carlo simulations 
using the QGSJET model as generator. 
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